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A precise modeling of Phoebe's rotation 
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ABSTRACT 

Aims. Although the rotation of some Saturn's satellites in spin-orbit has already been studied by several authors, this is not the 
case of the rotation of Phoebe, which has the particularity of being non resonant. The purpose of the paper is to determine for 
the first time and with precision its precession-nutation motion 

Methods. We adopt an Hamiltonian formalism of the motion of rotation of rigid celestial body set up by Kinoshita (1977) based 
on Andoyer variables and canonical equations. First we calculate Phoebe's obliquity at J2000,0 from available astronomical data 
as well as the gravitational perturbation due to Saturn on Phoebe rotational motion. Then we carry out a numerical integration 
and we compare our results for the precession rate and the nutation coefficients with pure analytical model. 
Results. Our results for Phoebe obliquity (23°95) and Phoebe precession rate (5580" .65/cy) are very close to the respective 
values for the Earth. Moreover the amplitudes of the nutations (26" peak to peak for the nutaton in longitude and 8" for the 
nutation in obliquity) are of the same order as the respective amplitudes for the Earth. We give complete tables of nutation, 
obtained from a FFT analysis starting from the numerical signals. We show that a pure analytical model of the nutation is not 
accurate due to the fact that Phoebe orbital elements e, M and Ls are far from having a simple linear behaviour. 
Conclusions. The precession and nutation of Phoebe have been calculated for the first time in this paper. We should keep on the 
study in the future by studying the additional gravitational effects of the Sun, of the large satellites as Titan, as well as Saturn 
dynamical ellipticity. 
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1. Introduction 

Understanding the rotation of gravitationally interacting 
bodies with an ever increasing accuracy is an important 
challenge of celestial mechanics. Many authors have mod- 
elled with very good accuracy the rotation of the Earth 
as Woolard (1953), Kinoshita (1977). The other planets of 
the Solar system followed logically (see for instance Yoder, 
1997; Rambaux, 2007). Today thanks to new and very 
precise observational astrometric techniques the dynami- 
cal study of celestial bodies extended to comets, asteroids 
and satellites (Meyer, 2008; Sinclair, 1977; Wisdom, 1984, 
1987). 

Nevertheless a few studies have been made on the ro- 
tational motion of Phoebe, the ninth satellite of Saturn, 
since its discovery in 1899 by W.H Pickering. In 1905, 
F.Ross etablished that this motion is retrograde and gave 
for the first time the orbital elements refering to the 
mean equinox and ecliptic of date. Jacobson (1998) de- 
termined new orbital elements thanks to Earth based as- 
trometric observations from 1904 to 1996. More recently 
Emelyanov(2007) elaborated ephemerides of Phoebe from 
1904 to 2027. Others studies have been made on the com- 



position and the inertial parameters of Phoebe (Aleshkina 
et al., 2010). Phoebe is the only non-synchronous satellite 
of Saturn with rather well known physical parameters, as 
the moments of inertia. As Phoebe is a dissymetric body 
(large dynamical flattening and large triaxiality) it is in- 
teresting to study its rotation. 

Here for the first time, we propose to determine the 
combined motion of precession and nutation of Phoebe 
considered as a rigid body After presenting the model 
(Sect(2), we describe the orbital motion of Phoebe by fit- 
ting the curves of the temporal variations of the orbital el- 
ements a, e, M and L s with polynomial functions. Thanks 
to a fast Fourier transform (FFT) , the large periodic vari- 
ations around the mean elements are analysed and show 
that the orbital motion of Phoebe is far from being keple- 
rian (Sect|5J). Then we determine the motion of precession 
and nutation of Phoebe by numerical integration using 
the value of the obliquity obtained in section @] Finally, 
after quantifying the validity of the developments done by 
Kinoshita (1977) for the Earth, when applied to Phoebe, 
an analytical model is constructed. 

All our results are compared to the results obtained 
by Kinoshita for the Earth. We show in section [5] that the 
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numerical approach is significantly more accurate to de- 
termine the precession and the nutation motion of Phoebe 
than the analytical model on the opposite of the Earth for 
which a very good agreement is reached. This is due to the 
large periodic departure from the keplerian motion of the 
orbit of Phoebe. We note however that it gives a fairly 
good approximation for the precession which is sensitive 
to the average value. 

2. Theoretical equation of the precession-nutation 
motion 

The present study of the rotation of Phoebe considered as 
a rigid body, is based on the theoritical framework set by 
Kinoshita (1977) who adopted an Hamiltonian formula- 
tion using Andoyer variables (1923) to evaluate the motion 
of the precession-nutation of the rigid Earth. Here a sim- 




Fig. 1. The precession nutation motion 

plificd model is described and only the variables and the 
equations used in this paper are detailed. The rotation of 
a rigid Phoebe involves three axes : the figure axis, which 
coincides with the axis of the largest moment of inertia, 
the angular momentum axis directed along (f, with G the 
amplitude of the angular momentum and a reference axis 
arbitrarly chosen as the axis perpendicular to the orbit 
of Phoebe at a given epoch t (see FigO}. The precession 
and nutation motions are respectively the linear and the 
quasi-periodic parts of the motion of the figure axis or of 
the angular momentum axis with respect to the reference 
axis. Assuming that the angle J between the angular mo- 
mentum axis and the figure axis is small as it is the case 
of the Earth (for which its value is less than 1"), only 
the motion of the angular momentum axis is considered 
here. The motion of this axis with respect to the refer- 
ence axis is described by the angles h,I (Andoyer, 1923) 
where / is the obliquity angle and h is characterizing the 
precession-nutation in longitude which corresponds to the 
angle between the reference point 7 t and the node Q. The 
reference point 7< is the intersection between the orbit and 
the equator of Phoebe at the date t , so-called " departure 
point" (Capitaine, 1986) and Q is the ascending node be- 
tween the plane normal to the angular momentum and the 
orbital plan. The Andoyer variables g, I in the Fig[T] are 



respectively the angle between the node Q and the node 
P and the angle between a meridian origin and the node 
P where P is itself the ascending node between the plane 
normal to the angular momentum axis and the equatorial 
plane. The proper rotation of Phoebe is described by the 
angle l + g = $. The Hamiltonian related to the rotational 
motion of Phoebe is: 



K = F a + E + E' + U 



(1) 



F Q is the Hamiltonian for the free rotational motion de- 
fined by : 



_ 1 sin 2 / 



^)(G 2 -£ 2 ) + ^. (2) 



where L is the component of the angular momentum axis 
along the figure axis and A, B, C are the principal mo- 
ment of inertia of Phoebe. E + E 1 is a component related 
to the motion of the orbit of Phoebe, which is caused by 
planetary perturbations and have been given in detail in 
Cottereau and Souchay (2009) when studying the rota- 
tion of Venus. U is the disturbing potential due to Saturn 
considered as a point mass and its disturbing potential is 
given by : 

u = U 1 + U 2 



GAP 

A-B 



2C-A-B 



P 2 (sin<5) 



P 2 2 (sin<5) cos 2a 



(3) 



^ GM'Mpa" r . 

n—3 
n 

— P™(sin<5) ■ (C nm cos ma + S nm sin ma)] (4) 

m— 1 

where G is the gravitational constant, M' is the mass of 
Saturn, r is the distance between its barycenter and the 
barycenter of Phoebe, a and S are respectively the longi- 
tude and latitude of Saturn (not to be confused with the 
usual equatorial coordinates), with respect to the mean 
equator of Phoebe and a meridian origin. The P™ are the 
classical Legendre functions given by: 

(-l) m (l - x 2 )^ d n+m (x 2 - l) n 



(5) 



2 n n\ d n+m x 
Notice that the perturbations due to the other planets, 
satellites and the Sun a priori of second order, will not be 
studied in this paper. The variation of the obliquity and 
the precession angle is given by : 



1 , 1 dK 



dl 
~dt 
dh 

~dl GsmI dl 



G sin/ dh 
1 dK 



cot / 



dK 



(6) 



To solve these equations, Kinoshita (1977) used the Hori's 
method. As the order of the disturbing function is : 

^ - M ' ^ 2r 2C-A- B) ^ 346 . 10 _ 8 (7) 



F M' + Mpto 



eyc- 



2C 



L. Cottereau, E. Aleshkina and J. Souchay: A precise modeling of Phoebe's rotation 



3 



the same method can be applied to Phoebe. Neglecting 
the very small contribution as the component E + E' and 
applying the first order of the Hori's method, this yields : 

A/ = U^llL I Utdt-cotl-^- flJidt] (8) 
G sin 1 oh J og J 



Ah=- 



1 







Gsinldl 



Uidt. 



(9) 



U\ given by ([3]) can be expressed as a function of the 
longitude A and the latitude (3 of Saturn with respect to 
the orbit of Phoebe at the date t using the transformations 
described by Kinoshita (1977) and based on the Jacobi 
polynomials: 



Ui = 



GM' 2C-A-B/ 1. 2r 
( --(3 cos 2 1-1) 



-- sin 2 Jcos2(A - h) 
4 

A — B 



-sin 2 Icos(2l + 2g) 



V -(1 + ecos/) 2 -cos2(A 

£ = ±1 



h — el — eg) 



(10) 



Finally the motions of precession-nutation in longitude 
and in obliquity are given starting from ([5]), ^ and JTI 
by 



Ah 



AI = 



1 d 

sin / d I 

1 d 
sin / dh 

d 



cot I- 



dg 



(Wi + W 2 )dt 



(Wi + W 2 )dt 



W 2 dt 



(11) 



(12) 



where 



GM' / 2C - A - B r 1 , , 
Wi = — ( - [--(3cos 2 /-l) 



Gr 3 



-- sin / cos 2 

4 



(X-h)]) 



(13) 



Wo 



GM' ( A-B 



Gr 3 



3 

-sin 2 /cos(2/ + 2g) 



+ Y / -(l + ecosir 

e=±l 

cos 2(A — h — el — eg) 



(14) 



Since W\ does not depend on g, this variable does not 
appear in the second part of the right hand side of equa- 
tion (|12p . Supposing that the components uii and ui 2 of 
the rotation of Phoebe are negligible with respect to the 



component W3 along the figure axis, as it is the case for 
the Earth, this yields : 



3GM' /2C~A-B r 1,, 



1 



sin 2 /cos2(A - tij] 



(15) 



W 2 



3GM' (A-B 



4G 

+ E J(l + "os/) 2 

e=±l 

cos 2(A — h — el — eg) 



2 sm2 I cos(2l + 2g) 



(16) 



where 



2C-A-B 



and 



A-B 



— cuiu iC are respectively the dynamical 
flattening and the triaxiality of Phoebe. As the rotation 
of Phoebe is retrograde the conventional notations of the 
precession and obliquity will be respectively : 



Atfj = Ah 
Ae = AL 



(17) 
(18) 



which are the opposite of the conventional ones for the 
Earth (Kinoshita, 1977, Souchay et al, 1999). 

3. Phoebe Ephemerides 

The ephemerides used in this paper for the Phoebe or- 
bital motion around Saturn have been constructed by 
Emelyanov (2007) where the orbital parameters and the 
coordinates (X, Y, Z) of Phoebe are given in the planeto- 
equatorial reference frame considered as inertial from the 
first observation date in 1904 to 2027. The reference point 
is one of the intersections between Saturn and Earth equa- 
torial planes at J2000.0. To apply the theoretical frame- 
work of Kinoshita (1977) the mean elements for the vari- 
ables a, e, M and L s are needed. In this section, these el- 
ements are calculated in the same way as Simon et al. 
(1994) for the planets. The differences betwen the true 
motion of Phoebe and a keplerian motion are also studied. 
Throughout the paper our study is restricted to the time 
span ranging from 2000.0 to the limit of the ephemeris 
2027.0. 

3.1. The semi-major axis 

The mean element for the semi-major axis for the Earth is 
a constant (Simon et al, 1994) and the relative variations 
around this mean value are of the order 10 ~ 5 . In compara- 
ison, for Phoebe, approximating a by a constant does not 
yield a good accuracy. FigJ2] shows the relative variation 
of a around its mean value estimated at a = 0.0864273 in 
9000 days time span. The variations consist clearly in peri- 
odic components with amplitude of the order of 10~ 3 . The 
leading oscillations of the signal are determined thanks to 
a fast Fourier Transform (FFT). Table Q] gives the lead- 
ing amplitudes and periods of the sinusoids characterizing 
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Fig. 2. Relative variation of a around the mean value 
a = 0.0864273 in 9000 days time span. The central curve 
represents the residual after substracting the sinusoidal 
terms of the Table [TJ 

Table 1. Leading amplitudes and periods characterizing 
the semi-major axis. 



M w 58.9233 + 0.656386- 1 - 0.000017- i 2 

+9.78065 • 10~ 9 • i 3 - 2.1541 • 10~ 12 • t 4 

+2.06592 • 10" 16 • t 5 - 7.2487 • 10" 21 • t 6 (20) 

L s 331.109 + 0.653712- 1- 3.64487- 10~ 7 • i 2 
+3.54186 • 10~ 10 • i 3 - 9.42668 • 10~ 14 • t 4 
+1.03032 • 10" 17 • i 5 - 4.02263 ■ 10~ 22 • t 6 (21) 

where t is counted in julian days, M and L s being ex- 
pressed in degrees. Fig|3] shows the variations of the ec- 
centricity of Phoebe for a 9000 days time span. An im- 
portant feature of Phoebe orbit is the large eccentricity 
as well as its very large relative variations (of the order 
of 10%). As we will see in section [6] this will be very im- 
portant in our study because of the large corresponding 
variations of (^) 3 that it causes for the orbit. 



Period amplitude amplitude 
sin cos 
days AU AU 



261.98 -2.319 1(T 4 2.609 10 -4 
498.93 -3.309 10 -5 -1.756 10 -4 
177.54 -5.715 10 -5 -9.971 10 -6 



the signal in FigJ5] where the central curve represents the 
residuals after substraction of these sinusoids. 

The presence of the relative large periodic components 
due to the additional perturbations on the orbit of Phoebe 
shows that the motion of this satellite is not close to a 
keplerian one, as it is the case for the planets of the Solar 
system. The effects of the attraction of the Sun, Jupiter 
and the other satellites of Saturn on Phoebe 's orbit entail 
an important departure from the keplerian motion. The 
individual study of the other orbital elements confirms 
this result. 

3.2. The mean elements of the orbital parameters 

In the following, the formula for the mean elements of 
Phoebe keplerian motion e, M and Ls are given, where 
M, L s are respectively the mean anomaly and mean lon- 
gitude with respect to the " departure point" ■ They were 
obtained by fitting the curves of the temporal variations of 
these elements given by Emelyanov (2007) with a polyno- 
mial expression at order 6 as it has been done for the plan- 
ets of the Solar system (Simon et al, 1994). Considering 
the limits of the ephemerides, truncating the polynomial 
functions at the order 6 looks as a sufficient approxima- 
tion. Thus : 

e w 0.168480 - 0.0000594-4 + 4.931508- 10" 8 ■ t 2 
-1.30162 ■ 10~ n ■ t 3 + 1.1744 ■ 10" 15 • t 4 
-5.2502 • 10~ 22 • i 5 - 3.0861 • 10" 24 ■ i 6 , (19) 



Fig. 3. Variation of the eccentricity of Phoebe for a 9000 
days time span. The central curve represents the polyno- 
mial function given by (1191) . 




Fig. 4. Residuals after substraction of the polynomials 
given by (|T9)l from the eccentricity of Phoebe for a 9000 
days time span. The central curve represents the residuals 
after substracting the sinusoidal terms of Table O 

Figures 21 [5] and |(5] represent respectively the residu- 
als after substraction of the polynomial functions given by 
(flTJl) . (|20p and (|2"l~T) from the eccentricity, the mean anomaly 
and the mean longitude of Phoebe. As already observed 
for the semi-major axis, the residuals have periodic com- 
ponents which reach an amplitude of 4° for the mean 
anomaly. 
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Fig. 5. Residuals after substraction of the polynomials 
given by (j2T))) from the mean anomaly of Phoebe M for 
a 9000 days time span. The central curve represents the 
residuals after substracting the sinusoidal terms of Table 
® 




of Phoebe is determined both analyticaly and numericaly. 
The differences between the results given by the two meth- 
ods due to the non keplerian motion of Phoebe will be 
significantly increased when compared with the equivalent 
study for the Earth (Kinoshita, 1977). All the Phoebe pa- 
rameters involved in the equations ([TTj) and ([T2"j) which 
enables one to calculate the precession and nutation are 
known, at the exception of the obliquity of Phoebe at 
J2000.0, which represents a fundamental initial condition. 
It is calculated in the following section. 

4. Obliquity of Phoebe at J2000.0 

Using a similar method as in Cottereau and Souchay 
(2009), the obliquity is given by : 

cos I Q = P^ P -P^ (22) 

where P op and P p are respectively the unit vectors directed 
toward the orbital pole and the pole of rotation of Phoebe 
at to, defined by : 

(sin iq sin fio \ 
— sinzocos^o (23) 
cos« / 



Fig. 6. Residuals after substraction of the polynomials 
given by (f2"Tj) from the mean longitude of Phoebe L s for 
a 9000 days time span. The central curve represents the 
residuals after substracting the sinusoidal terms of Table 

m 



3.3. Frequencies analysis 

Thanks to a FFT analysis the leading ampitudes and 
the periods of the sinusoids characterizing the signals in 
Figs|H[S]and|n]are determined and given in Tabled] 



Table 2. Leading amplitudes and periods characterizing 
the orbital elements e, M and L s . 



orbital element 


days 


amplitude 


amplitude 




501.078 
3505. !>9 
550.62 


7.1200 10" 4 
1.4949 10~ 3 
1.260 10" :i 


5.776 10 _:i 
8.19196 10 -4 
-8.9235 1<T 4 


Mean anomaly 
M 


498. 88 
553. ;« 
177.50 


2.4850 
0.5212 
1.2010" 2 


-0.488 
0.2515 
-0.2338 


Mean longitude 


262.02 
499.25 
554.2945 


-0.2959 
0.3470 
6.727210" 2 


-0.2669 
-5.211610" 2 
7.101110" 4 



I COS Sp cos 



n 



sm a" cos 



V 



sin 5° 



(24) 



/ 



io, S7o are the inclination and the longitude of the as- 
cending node of Phoebe in the Saturn equatorial reference 
frame at J2000.0 given by the ephemerides of Emelyanov 
(2007) and a®, 5® are the geo-equatorial coordinates of the 
north pole of Phoebe given by Seidelmann et al. (2007). 
The numerical values of these variables are given in ap- 
pendix. Applying two rotations to P p , the two units vec- 
tors are expressed in the Saturn equatorial reference frame 
and the obliquity at J2000.0 can be determined directly 
from (I2"2"j) . 

After computation and taking into account the Phoebe 
retrograd rotation the numerical value of the obliquity at 
J2000.0 is : 

I Q = 23°. 95 (25) 

Notice that this obliquity is close to the Earth 
one (23° .43). Knowing the obliquity at J2000.0, the 
precession- nutation of Phoebe can be determined with the 
help of (fTTj) and (|12p in the following sections. 



The central curves in Figs 01 [5] and [5] represent the 
residuals after substraction of the leading sinusoidal com- 
ponents of the mean elements. We want to highlight here 
the fact that the motion of Phoebe is far from being a 
quasi keplerian motion as it is the case for the planets of 
the Solar system. In the following the precession-nutation 



5. Numerical results of the precession-nutation of 
Phoebe 

Using the ephemeris of Emelyanov (2007), the precession 
and the nutation in longitude and in obliquity are de- 
termined by numerical integration with an extrapolation- 
algorithm based on the explicit midpoint rule. We remind 
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here that the precession and the nutation in longitude and 
in obliquity are denned as the variation in space of the an- 
gular momentum axis with respect to the reference axis 
due to Saturn's gravitational perturbation. They corre- 
spond respectively to the linear part of h and the quasi- 
periodic part of h and /. Before integration, the planeto- 
equatorial coordinates (X, Y, Z) taken from ephemerides 
are converted into rectangular coordinates with respect to 
the orbital plane, which at their turn are converted into 
the spherical coordinates used in (fTTj) and (fT2l . The values 
of the dynamical flattening and the triaxiality of Phoebe 
are taken in Aleshkina et al. (2010). Figs J7] and |S] show 
respectively Phoebe precession-nutation in longitude and 
Phoebe nutation alone for a 9000 days time span (24.6 y). 
FigJS] shows the nutation in obliquity for a 2000 days time 
span (5.47 y). Notice that 9000 days correspond to the 
limit of the ephemerides. The 2000 days time span of the 
nutation in obliquity is chosen in order to see clearly the 
leading oscillations. The precession and the leading oscil- 
lations of the nutation in longitude and in obliquity are 
determined thanks to a linear regression and a fast Fourier 
Transform (FFT). Our value obtained for the precession 
of Phoebe is : 

* = 5580.65* (26) 

where t is counted in Julian centuries. Notice that this 
value for the precession in longitude is very close to the 
precession for the Earth ($ = 5000" .3cy). So the effect of 
the tidal torque of Saturn on the precession of Phoebe is 
roughly the same as the combined gravitational effect of 
the Moon and Sun on the precession of the Earth. The 
physical dissymmetry (large values of the dynamical flat- 
tening and of the triaxiality) of Phoebe which directly 
increase the amplitude of the precession are compensated 
by its slow revolution and fast rotation (see eq (|15p ) . 




Fig. 7. The precession and the nutation of Phoebe in lon- 
gitude for a 9000 days time span. 

The nutation in longitude is given by : 

A* = V" $j sinfat) + $ ■ cos(^i) + ®it sxD.(i/it) 

+9-tcos(^t) (27) 

where $i, <£>-, 6i, 9- and Vi are respectively the lead- 
ing amplitudes and the frequencies of the sinusoidal terms 




2000 4000 6000 8000 



Fig. 8. The nutation of Phoebe in longitude for a 9000 
days time span. The central curve represents the residual 
after substracting the sinusoidal terms of the Table |3] 
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Fig. 9. The nutation of Phoebe in obliquity for a 2000 
days times span. The central curve represents the residual 
after substracting the sinusoidal terms of the Table 2] 

and can be found in Table [21 Similarly, the nutation in 
obliquity is given by : 

AI = 2J Ti sin(/i,t) + r • cos(^t) -I- ILt sin(/i,t) 

+n^cos(/Xit) (28) 

where r,, T^, IL,, IT^ and are given in Table [4] The corre- 
sponding arguments as a combination of Phoebe mean lon- 
gitude and mean anomaly are determined empirically and 
given when they have been clearly identified. The residu- 
als after substraction of (|2"T1) and (J2HJ) to the results of the 
numerical integration are the flat curves in Figsj8]and[9l 

Table 3. A^ = Ah: nutation coefficients in longitude of 
Phoebe 



Argument 




arc second 


** 

arc second 


isin 

e 

uc second (10 -6 ) 


;cos 
& 

arc secondUO- 15 ) 


2L S 


275.1353 


-3.411 


5.8789 


20.7222 


-5.3705 


2L S - M 


550.721 


3.5260 


4.29957 


48.644 


-44.6566 


2L S + M 


1S3.44 


-2.4496 


0.2363 


11.748 


21.716 




5564.9 


-0.2884 


0.52829 


-63.3917 


31.786 


2L S + 2M 


137.6 


-0.482874 


- 0.5341 


-5.0143 


10.129 




497.34 


-0.5300 


0.14200 


18.758 


8.77425 


2L S + Ml 


110.09 


0.0348 


-0.1723 


-4.1343 


-2.5833 


2L S - 5M 


177.17 


0.0413 


-0.2188 


-1.4336 


7.4353 


2L S - QM 


134.04 


0.092656 


-0.00384 


-3.1435 


0.243 




3220.0 


0.1203 


0.1494 


-1.53677 


-30.596 




260.00 


-0.1160 


-0.0257 


13.69 


2.3207 
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Table 4. Ae = Al : nutation coefficients in obliquity of 
Phoebe 



Argument Period 



II' 



second 



p second 



md(l(r' 



2L S 
2Ls + M 
2Ls - M 
ILs + 2M 
2L„ + 3M 
2Ls - 5M 
2L S - 6M 



275.12 
183.45 
549.90 
137.60 
110.09 
177.15 
134.04 

8166 
614 

91.75 
291.67 
107.80 
203.03 



-2.266 
-0.0937 
0.435 
0.2326 
0.0819 
0.0910 
0.01659 
-0.0089438 
-0.025055 
0.0078745 
0.010783 
-0.007745 
0.006117 



-1.4174 
-1.03465 
-0.1703 
-0.2053 
0.0164 
0.01957 
0.0456 
-0.1403 
-0.02073 
0.01760 
-0.0228555 
0.0123306 
-0.014634 



14786 
-9.01515 
0.92524 
-4.2893 
-0.2074 
-3.3764 

0.1052 
-5.4749 

1.1722 
0.72270 
-1.17099 
0.42834 
-4.14544 



1.2152 
4.1209 
3.7466 
-2.4026 
-1.5682 
-0.7907 
-1.4259 
19.3857 
-0.7703 
-0.1035(1 
1.0911 
-0.0620 
3.24460 



The nutation in obliquity is dominated by two fre- 
quencies associated with the arguments 2L S and 2L S + M 
with respective periods 275.13 d and 183.45 d. In addi- 
tion to these two frequencies, the nutation in longitude 
also presents another leading component with argument 
2L S — M. The components which are not identified in 
Tables and H as those with period 497 d and 260 d, 
already pointed out in Tables Q] and [2] (see SectJ3]). They 
probably come from the non Keplerian motion of Phoebe. 
As in the case for the precession, the nutations in longitude 
and in obliquity of Phoebe, which show respectively peak 
to peak variations of 26" or 8" are of the same order as the 
nutation of the Earth (respectively with a factor | and i). 
As the value of the triaxiality is large (Aleshkina et al., 
2010), it is interesting to evaluate its effect on the nutation 
of Phoebe. This effect is determined by the equations (TTTT) 
and (fl"2"l) taking only into account the component depend- 
ing on in equations (HU) and (HU). Figs [TU1 and [TTl 
show respectively the nutation in longitude and in obliq- 
uity depending on the triaxiality for a 10 days times span 
beginning at 410 days. This short time span is chosen to 
see clearly the high frequency oscillation to the rotation 
of Phoebe with period roughly half the proper rotation 
of Phoebe (9 hrs) and the higher amplitude of the whole 
signal. Thus we remark that despite the large triaxiality 
of Phoebe, its effect on the nutation in longitude and in 
obliquity, at the level of a few 10 -3 arcseconds amplitude, 
is very small . As explained already for the precession, the 
effect of the disymmetry of Phoebe, characterized by a 
rather large value of its triaxiality, is compensated by its 
very fast rotation, which strongly decreases the amplitude 
of the signal after integration. 

The goal now is to determine an analytical model to 
describe the precession-nutation of Phoebe and to com- 
pare it with our numerical results. In the following the 
effect of the triaxiality on the nutation will be ignored 
and the in (0 and CO} are replaced by ^f-{f) 3 . 

6. Construction of an analytical models for the 
precession-nutation of Phoebe 

Many authors have elaborated an analytical model of the 
rotation of the Earth in order to understand the causes 
of the observed motion and to have the ability to predict 




















Fig. 10. Effect of the triaxiality on the nutation in longi- 
tude for a 10 days time span. 




Fig. 11. Effect of the triaxiality on the nutation in obliq- 
uity for a 10 days time span. 



it (Kinoshita, 1977, Bretagnon et al, 1997, Roosbeek and 
Dchant, 1997, Souchay et al., 1999). In this approach, it 
is interesting to build an analytical model based on the 
model of Kinoshita for the Earth to determine the preces- 
sion and the nutation of Phoebe. 

To solve analytically the (fTTj) and (|12"j) . Kinoshita 
(1977) developed \{^f and (^) 3 cos2A as a function of 
time, through the variables M and L s taking the eccen- 
tricity as a small parameter : 

I(^) 3 = /(M,e) + ( e 3 ) 
2 r 



(-) 3 cos2A = g(M, L s ,e) + o(e 3 ) 
r 



(29) 



The detailed demonstration is presented in Kinoshita 
(1977). As the eccentricity of the Earth is small, trun- 
cating the development at the 3rd order is a sufficient ap- 
proximation. On the opposite the eccentricity of Phoebe 
is higher (by one order) than the eccentricity of the Earth 
(0.18 instead 0.0167). Moreover, as we have seen before, 
Phoebe orbital elements are undergoing large variations. 
Then the relative error made by replacing 5(f) 3 and 
(^) 3 cos2A by the developments such as (|2T)1) need to be 
quantified. To that purpose, \{^) 3 and (^) 3 cos2A are 
developed as for the Earth but with the real values of 
e,M,Ls taken in the ephemerides of Emelyanov (2007), 
instead of a linear value (for e, M ans L s ) for the Earth. 
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6.1. Test for the developments of ^(^) 3 and 
(^) 3 cos2A. 

In this section, the developments of an d (;:) 3 cos2A 

done by Kinoshita (1977) in the case of the Earth and 
by Cottereau and Souchay (2009) in the case of Venus 
are tested for Phoebe. We used the same methods as in 
Cottereau and Souchay (2009) to determine the develop- 
ment of T^if) 3 - As the eccentricity of Phoebe is large, the 
development at the 3rd order as made by Cottereau and 
Souchay (2009) for Venus is insufficient as it yields rela- 
tive errors to the order of 4.10~ 3 . This is due to the slow 
convergence of the development which must be carried out 
up to 6th order in our case, if we want to reach a relative 
10 -4 accuracy. FigfT2lshows the discrepancies between the 
numerical and the semi analytical results using this devel- 
opments both to the 3rd and 6th order. This shows that 
high orders developments are needed when taking into ac- 
count highly eccentric orbits as Phoebe's one. 




Fig. 12. Ai(^) 3 between the numerical and the semi an- 
alytical results using the developments to 3rd (blue curve) 
and 6th order (red curve). 



3 cos2(A - h) = (1 - \? + H e 4 + ^e 6 )cos2L s 



4 81 



e 5 cos(2L s - 5M) 



+ (e 4 ^ + ^e 6 )cos(2L s -4M) 
+ (e 3 ^ + ^e 5 )cos(2L s -3M) 

+ (-^ + e 3 l-A^)cos(2L s -M) 

.7 2 3 123 489 5 . 
+ ( 2 C +e V + 128 e)C ° S(2is 

+ (^e 2 + e 4 ^ + ^e 6 ) cos(2L s + 2M) 



M) 



.845 3 

+ ( e 3 

V 48 



,32525, 
768 ' 



,533 4 fi 13827, 
+ ^!l! e 5 cos (2L s + 5M) 



cos(2L s + 3M) 
cos(2L s + 4M) 



3840 
73369 

720 



e 6 cos(2i s + 6M) + 0(e 7 ) 



(31) 




For the developments of (^) 3 cos 2 A, the same conclu- 
sion is obtained. FigJT3] shows the discrepancies between 
the numerical and the semi analytical results using the 
developments of - 3 cos 2 A, to 3rd and 6th order. The dif- 
ference is descreased from 0".3 to 0".01 (peak to peak). 
The developments are given by : 



2 V 



1 3 2 
= 2 + 4 6 + 

.3 27 3 
+ ( 2 e+ 16 6 



15 
16' 



261 
128 



35 6 

e b 

32 

e 5 )cosM 



-(7 e2 + 7 e4 + ^ e6 ) cos2M 
v 4 4 64 ; 



3167 
+ ~320 



393 
256* 
129 

160 



-e b )cos3M 



e b )cos4Af + 



1773 
256~ 



-e° cos5M 



e 6 cos6Af + 0(e 7 ) 



(30) 



Fig. 13. A(^) 3 cos 2A between the numerical and the semi 
analytical results using the developments to 3rd (blue 
curve) and 6th order (red curve). 



6.2. Analytical model of Phoebe 

To express and (^) 3 cos2A by purely analytical 

functions, Kinoshita replaced e,M,Ls by their linear 
mean values. For Venus, the same method is used in 
Cottereau and Souchay (2009), who replaced e, M, Ls by 
their linear mean values given by Simon et al (1994). Using 
the mean elements of the section ([2]) and the developments 
given by (|30l) and (|3ip . an analytical model is obtained. As 
explained in Kinoshita and Souchay (1990) and Cottereau 
et al. (2010) for the Earth and for Venus, the difference 
between the nutation computed by analytical tables and 
the numerical integration is small (at the order of 10~ 5 ) 
and caused by the indirect planetary effects on the planet. 
For Phoebe, this difference is also caused by indirect per- 
turbations but cannot be ignored as it will be presented 
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in the following. Replacing the semi-major axis in (1111) 
and (|T2|) by a (Sect [3]), Ah and Al depend directly on the 
scaling factor K$ ■ 

K s = —^-^r- - = -11674" .4857 Julian.cy 

(32) 

where 2C 2 ^ B is the dynamical flattening of Phoebe, w 
the angular velocity of Phoebe and n the mean motion of 
Phoebe given by : 

n 2 a = GM' (33) 

The numerical value is deduced directly from the values 
of the A,B,C and u taken respectively in Aleshkina et 
al. (2010) and Bauer et al (2004). Notice that this scaling 
factor is very close to the Earth one (11036". 6 by Julian 
cy). As seen in the section ([2]) the precession of Phoebe 
corresponds to the polynomial part of (lll[) and can be 
directly determined by : 

V, =-(^cos/|(i + ^ 2 + ^e 4 + |e 6 )^) 

= 5541.18* (34) 

where e is replaced by the linear part of (|19[) . The discrep- 
ancies between the precession as computed by analytical 
integration and that given by (|2"rJ|) are of the order of 1% 
and show that the analytical model describes with a good 
accuracy the precession of Phoebe. Replacing respectively 
e, M, Ls in the developments ((30]) and ([31]) by the linear 
part of (|T9" ]) .([2"0" f and ([2T|). the nutation of Phoebe is deter- 
mined by analytical integration. Figures [JJ] and Qj)] show 
the differences between the nutation in longitude and in 
obliquity of the angular momentum axis from a purely 
analytical formulation and numerical integration on 9000 
days time span. These differences are compared to the 
residuals obtained in section [5] after substraction ([2~T1) and 
(|2"8"|) determined by FFT to the numerical signal (central 
curves in the Figures HH and [15]) . 




Fig. 14. Difference between the nutation in longitude of 
the analytical and numerical integration for 9000 days 
time span. The flat curve represents the residual after sub- 
straction (|2"T|) to the nutation signal in longitude obtained 
by numerical integration. 

Figures [T4l and [TBI show that the semi- analytical func- 
tions characterizing the nutation of Phoebe respectively 




Fig. 15. Difference between the nutation in obliquity of 
the analytical and numerical integration for 9000 days 
time span. The flat curve represents the residual after sub- 
straction (j2"8")) to the nutation signal in obliquity obtained 
by numerical integration. 

in longitude and in obliquity determined thanks to a fast 
Fourier Transform are more accurate than these obtained 
by a pure analytical integration (with a linear values of 
e,M,L s ). The pure analytical model restricts the expres- 
sions of the functions which describe the orbital elements. 
For the Earth (Kinoshita and Souchay, 1990) and for 
Venus (Cottereau et al. 2010), this model is very satis- 
fying because their motion are close to a keplerian one 
and consequently the orbital parameters can be replaced 
by linear functions. As the motion of Phoebe is highly dis- 
turbed in comparison to the planets of the solar system, 
the curves of the temporal variations of the orbital ele- 
ments (Emelyanov, 2007) cannot be fitted by simple func- 
tions as shown by the large variations around the mean 
elements in Sect [31 Thus the analytical model does not de- 
scribe the periodic components due to the indirect effect 
of the Sun or of the other satellites of Saturn. In contrast 
these effects are directly taken into account when using 
the functions given by ([2~T|) and (|28p because the FFT ap- 
proach followed in Sect[5Jis better suited to describe quasi- 
periodic motion. However the precession is less sensitive 
to these periodic variations, which are averaged during the 
calculation. This explains the similarity of the precession 
rate given by the two models. To conclude, the purely ana- 
lytical model set by Kinoshita (1977) for the Earth gives a 
good first approximation of the nutation of Phoebe but a 
more sophisticated analytical model is needed if one wants 
a precision down to the 10~ 5 level. We precise here that 
this accuracy is internal to our model and is independant 
of the uncertainties in the observations. This shows that 
we are limited only by the uncertainties on the input data 

7. Conclusion 

The purpose of this paper was to calculate for the first 
time the combined motion of precession and nutation of 
the Saturn satellite Phoebe. 

One of the important steps in our study was to describe 
the orbital motion of Phoebe by fitting the curves of the 
temporal variations of the orbital elements a, e, M and L s 
with polynomial functions. As shown in (3[ the periodic 
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variations of the semi-major axis around its mean value 
with a relative amplitude of the order of 10 -3 are clearly 
due to the perturbing effects of other celestial bodies like 
the Sun. The orbital motion is far from being a keple- 
rian one, as shown by the large polynomial expressions of 
e, M and L s , and by the large sinusoidal amplitudes char- 
acterizing the residuals after substraction of these poly- 
nomials. Applying the theoritical framework already used 
by Kinoshita (1977) for the Earth, the precession and the 
nutation motion of Phoebe are determined both analyti- 
cally and from numerical integration of the equations of 
motion. We found that the precession-nutation motion 
of Phoebe undergoing the gravitational perturbation of 
Saturn is quite similar to that the Earth undergoing the 
gravitational effect of both the Moon and the Sun. Thus 
our value for the precession of Phoebe, that is to say 
5580". 65 cy, is very close to the corresponding value for 
the Earth (5081" /cy) and the nutation in longitude and 
in obliquity of Phoebe with peak to peak variations of 26" 
and 8" are of the same order of amplitude as the nuta- 
tion of the Earth (respectively 36" and 18" peak to peak). 
Moreover Phoebe obliquity (23°. 95) is roughly the same 
as the Earth's one (23°. 43). Notice that the physical dis- 
symmetry (large value of the dynamical flattening and of 
the triaxiality) and the large eccentricity of Phoebe which 
direclty increase the amplitude of the precession and the 
nutation is compensated by its slow revolution and fast 
rotation. 

In Cottereau and Souchay (2009), we have demon- 
strated that the effects of the large triaxiality of Venus 
on its nutation are of the same order of magnitude than 
the effects of the dynamical flattening. By contrast, the ef- 
fects of the large triaxiality of Phoebe on its nutation are 
compensated by its fast rotation and decrease the ampli- 
tude which become negligible compared to the dynamical 
flattening. 

We also investigated the possibility to construct an- 
alytical tables of Phoebe nutation, as was done for the 
Earth. As Phoebe has a large eccentricity, the analytical 
developments at the 3th order given by Kinoshita (1977) 
for the Earth are not valid. We must carry out the devel- 
opments up to 6th order to reach a relative 10 -4 accuracy. 

Thus although the amplitude of the nutation motion 
is close to the Earth one, we demonstrated that the ana- 
lytical model used by Kinoshita for the Earth does not 
describe the nutation motion of a disturbed body like 
Phoebe with the same accuracy. This analytical model 
does not take into account the large perturbing effects of 
the celestial bodies on the orbit of the satellite. 

To describe the nutation motion of Phoebe a FFT ap- 
proach is better than a pure analytical integration done 
with linear expressions for e,M and L$ as shown Figs ll4l 
and [T5] The FFT analysis is better fitted to describe the 
periodic variations characterizing the nutation signals of 
Phoebe. These periodic variations are averaged during the 
calculation of the precession which decreases the discrep- 
ancies between the two models. 



To conclude we have shown that the analytical model 
set by Kinoshita (1977) gives a good first approximation of 
the precession-nutation of Phoebe but further analytical 
developments are needed to reach the same accuracy than 
for the terrestrial planets. 

We think that this work can be starting point for fur- 
ther studies such as the elaboration of another very pre- 
cise analytical model of the rotation of Phoebe by taking 
into account effects ignored in this paper, as the direct 
effects of the Sun, of Titan and of Saturn dynamical flat- 
tening. Such a model is required to develop the long term 
ephcmcrides of Phoebe's rotation, which should require 
long term orbital ephemerides, not still available. 
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Table 5. Numericals values used in this paper 



Phoebe 


Eccentricity (,12000.0) 


0-1648 


Inclination (.12000.0) 


151°. 64 


Ascending node (.12000.0) 


54°. 317 


Obliquity 


23° .95 


Period of rotation 


0.386396 d 


Triaxiality : 


-0.011125 


Dyn. flattening : ' 2C ~ 2 £~ B 


0.06465 


geo-equatorial coordinates of the pole of 


Phoebe l*£ = 356°. 90 and S° = 77°. 80 



siru-equat urial cuunlinat es ui I he pole ol Saturn aso — 40° .589 and Ssn = 83"' .537 



